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Abstract 


The present investigation involves ultrasonic tomographic non-destructive evaluation of 
composite specimens implanted with artificial inclusions of polystyrene in perspex and 
evaluation of effectiveness of different features of ultrasonic signal in detecting different 
types of flaws. 

Four different specimens are used which are made of perspex and the inclusions axe 
of polystyrene. The wedge is used only for the normal incidence of the ultrasonic ray 
so that the transmission through the specimen is maximum. A through-transmission 
technique is used for data collection. The test section is immersed in the water which act 
as a couplant. After collecting the data for the required number of rays and views and 
applying correction is given as the input to the Convolution Back Projection algorithm. 
The data is reconstructed with different class of filters such as cosine, sine, hamming and 
Ramachandnm-Lakhminarayan. 

The output of the reconstruction algorithm after pseudo colouring is an imago which 
show dearly the different acoustic impedences zones. So flaws can easily be detected. 

The results are quite encom aging. The present experimental set-up, after some mod- 
ification, can he used for detailed flaw detection and sizing in 3-Dimensional ultrasonic 
tomography. 
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N omenclat ure 


/ function reproMuit ing the spatial distribution of any 
physical quantity 
x, y cartesian coordinates 
7\ d polar coordinates 

p line integral data (projection data) 

N energy of the acoustic waave after travelling a distance dl 
Ao initial energy of the acoustic wave 
p slowness coefficient 
A Fourier transform 
R Fourier frequency 
R c cut-off frequency 
XV window function 

H f/# (0) second order derivative of window function 
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Chapter 1 


Introduction 


The technique of eomputrised tomography (CT) has established itself as a leading tool 
in diagonostic radiology over the past twenty years and it is gaining importance in the 
nondestructive evaluation area in a variety of situations. In conventional projection ra- 
diography, the result is obtained just in the form of a shadowgraph requiring subtle 
interpretation. CT, however, uses tomographic reconstruction to combine the informa- 
tion from multiple projections. Thus, the CT information is displayed quantitively as a 
reconstructed slice plant* of the object and it can provide accurate qualitative measures 
of material characteristics in small volume elements of the component. 

It is noted that tomographic extensions of classical measurement techniques incorpo- 
rate multiple projections, which increases the data collection time. The results, however, 
are much better if tomographic reconstruction is used. Tomography, if perfected, can 
also provide a cross validation of other experimental techniques. 
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1.1 Image Reconstruction from Projections 

The range of applicability of image reconstruction is very wide. At one end data from 
electron microscope are used to reconstruct the molecular structure of bacteriophages, 
while at the other end, data collected by the rockets sent outside the earth’s atmosphere 
are used to reconstruct the X-ray structure of supernova remnants. Of all the applica- 
tions, probably the greatest effect on the world at the largest has been in the area of 
diagnostic medicine. 

Image reconstruction can be very aptly defined in words of Herman(1980) as: 

“ Image reconstruction from projections in the process of producing an image of a 
two-dimensional distribution (usually of some physical property) from estimates of its 
line integrals along a finite number of lines of known locations”. 

1.2 An Overview of the process of CT 

The aim of computerized tomography (CT) is to obtain information regarding the nature 
of the material occupying exact positions inside the object, there are other modes of tra- 
ditional tomography, e.g. line tomography and transaxial tomography. In tradition form 
of tomography, objects which are out of focal plane are visible on the image, although 
in a blurred form. In CT, the images of cross-section are not influenced by the objects 
outside those sections. 


1.3 Composite Materials 

Demand imposed by modern advanced technologies on materials have become so diverse 
and severe that they often cannot be met by single component conventional materials. It 
is necessary to combine two or more materials to make a composite to meet the demand. 
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Constituents not only contribute their respective share, but their combined action tran- 
scends the sum of the individual properties and provides new performance, unattainable 
by constituents acting alone. Composites can be formally defined as materials having 
two or more chemically distinct materials or phases with clear interfaces on a macroscale 
among them. 

The fibre reinforced materials are probably the most important class of composite 
materials which are used widely in structural applications particularly when weight saving 
is crucial. These materials offer high strength-to- weight and modulus-to- weight ratios 
apart from their properties of excellent corrosion resistance, formability and mechanical 
properties which can be controlled to significant degree. 

1.4 Nondestructive Testing 

Since man and their machines are fallible, there is a need to inspect the components to 
assess the material. A wide variety of test schemes exist, some destructive and some non- 
destructive. The practical benefits of nondestructive testing are obvious by its name, as 
long as the results tire reliable and the inspection is cost-effective. In general, the various 
nondestructive evaluation (NDE) techniques can be placed into two categories: active 
and passive. The active techniques are those where something is introduced into or onto 
the specimen and a response is expected if a defect is present whereas passive techniques 
are those that monitor or observe the item during either a typical load environment or 
a proof cycle and attempt to determine the presence of a defect throgh some reaction 
of the specimen. Magnetics, ultrasonics, radiography fall in the category of active NDE 
and acoustic emission, noise analysis, leak testing etc. are the examples of passive NDE. 
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1.5 Ultrasonic Testing 

Ultrasonic inspection involves introducing a low energy, high frequency stress pulse into 
the material under inspection and examining the subsequent propagation of this energy. 
Defects such as cracks and inclusions act as a source of wave scattering throgh reflection, 
refraction, diffraction and mode conversion. It will be quite relevent at this stage to have 
a brief discussion on different types of waves. 

A wave is a transient phenomenon that depends upon time and position in space and 
which transports energy through space. The effect of a wave due to interaction with mat- 
ter may be measured by certain physical attributes such as displacement, strain, stress, 
density, electric field, magnetic field etc. In case of mechanical waves, it is convenient 
to classify them according to the geometry of the medium through which the wave is 
pi upag.it ing Thus, bulk waves are those that propagates through a medium so large in 
size compared to the wavelength that, the medium may be cosidered to be infinite. Bulk 
waves are generally plane or spherical waves depending on the shape of the wave front. 

Bulk waves in solids are normally described as being longitudinal (or congressional 
or dilatational) as the particle displacement caused at the points are parallel to the 
direction in which the wave is propagating. Sound waves in atmosphere is an example 
of this. A transverse (or shear) wave is one which produce particle displacements that 
are pvtpemiii ular to the direction of the wave propagation. Longitudinal waves typically 
travel at a faster speed than transverse waves. 

Surface or Rayleigh waves are observed in the vicinity of the boundary when bulk 
waves propagate parallel to the boundary in a semi-infinite solid medium. Its velocity 
is less than that of transverse waves in the same medium and it attenuates rapidly on 
entering the medium. The particle displacement produced by the passing of this wave is 
such that the path described by the particle over one complete cycle of the wave is an 
ellipse.These waves can be suitably used to detect surface flaws. 
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The interaction of bulk waves with geometrical constraints (in the form of a stress 
free parallel boundary, e.g., a plate) produces a very complicated structure of waves 
referred to as plate waves or Lamb waves. Lamb waves are of two types, symmetric 
and antisymmetric. In symmetric type, the particle displacement fields are symmetric 
with respect to the plate midplane and quite analogous to longitudinal modes. The lamb 
waves, in general, are composed of two components; standing weaves, that exist between 
the top and bottom surfaces of the plate and travelling waves, that propagate parallel to 
the midplane. 

For the present concern, the stress waves are generally assumed to be in the ultrasonic 
range, i.e., having frequencies greater than 20KHz. When an ultrasonic wave reaches an 
interface or a discontinuity, a portion of the energy is reflected and the other portion 
is transmitted. The amount of the reflected energy depends on the relative acoustic 
impedances of the media. Thus, the energy transmitted through the material is re- 
duced due to attenuations in the sample. Subsequently, these variations are reflected as 
the change in the relevant ultrasonic feature and serve as means for locating defects or 
discontinuities in the path of the wave. 

In practice, elastic waves are generated through the excitation of a peizocleetric device 
(in the probe) by electric pulse generator and the mechanical energy is then fed into the 
test specimen using a coupling medium (couplant) such as water, oil, grease etc. 

Depending on the mode of recieving of the waveform, testing are performed, in gen- 
eral, either in pulse-echo or in through-transmission mode. In pulse-echo mode of scan- 
ning, a high voltage spike sets the probe into oscillation and it becomes an open non- 
conducting electrical circuit and the reciever circuit waits for the return signal to strike 
the probe. In the through-transmission arrangement, the sending and the receiving end 
probes have different electronic circuits. Through-transmission is effective for detecting 
delamination , voids etc. In this case, however, both sides of the specimen should be 
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accessible. Pulse-echo technique can be employed if the part is not very thin, so that 
the initial pulse and the reflected pulse from the back surface do not overlap, nor the 
specimen is too thick, causing excessive attenuation by high signal loss even for defect 

free results. 

1*5.1 Presentation of ultrasonic data 

Hist orically, investigators have used the terminology A, B, C-scan to describe the type of 
p:« imitation of the data taken during ultrasonic study of specimens. An A-scan indicates 
a variation of signal amplitude at a point with time, where the ordinate of the oscilloscope 
trace* is proportional to the amplitude and abscissa to the time elapsed. A B-scan is the 
p: ‘iitation of the data in the form of a depth profile versus position along a cross section. 

A C-scan is the presentation of variation of any feature or feature set as a function of 
position of tin' transducer in two dimensionals relative to the surface of the examined 
•c t \ * Linen. During automated scanning, the transducer is generally translated along two 
mutually perpendicular axes in raster fashion using stepper motors. Essentially feature 
or feature set, used to construct the C-scan, is extracted from the recieved ultrasonic 
waveform obtained in the A-scan mode. Another recent way of representing ultrasonic 
data is by using tomographic technique which reconstructs a slice of an object from 
multiple projections. Thus, tomography, if perfected, can provide better results 

1.6 Reconstruction Algorithms 

Th<* simplifying assumptions in setting up the theory for reconstruction algorithms are 

1 . Slices are infinitely thin 

2. For any perticular source and detector pair position, wave travel in the same 
straight Jims (which lies in the infinitely thin slice). 
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A consequent ^ irs * Gumption is that the distinction between the voxel and 

the pixel disappears. Indeed, since the slice is infinitely thin, it can be thought of as 
a picture whose g* av I |>v el or the colour level at any point (x,y) is proportional to the 
relative slowness Todficiont /r(x,y) at that point. This is the reason why the theory be- 
hind reconst i net ion algorithm is referred to as “image reconstruction from projections”. 
Tomographic algorithms can he classified into two categories: 

1. Transform met heals 

2. Series expansion methods 

When the set of data available is large and complete, transform methods are generally 
used. Whereas if the* data set is not complete, series expansion methods are preferred. 
Transform methods an* riuert methods while series expansion methods can be direct as 
well as iterative, The scope of tin/ transform algorithm is limited to the reconstruction of 
a two-dimensional (2 D) image from integrals along a line in the plane, but they are useful 
in many application involving three-dimensional (3-D) objects becuase such object can 
he <1* comput'd into a stack of thin 2-D slices. If line integral data are available for each 
of these individual slices independent of others, then the 3-D reconstruction problem may 
he d* < omptiM-d into repeated 2 D reconstruction of serial cross-sections of the object. 


1.7 Ultrasonic Tomography 

Although CT was conceived originally for X-rays, in recent times, CT has been im- 
plemented successfully for other energy sources such as acoustic, laser, microwaves etc. 
Among these, acoustic tru:.».;*iaph> (especially using ultrasonic waves (UCT)), has spe- 
cial importance in apple ,»ti< m to NDE. Like X-ray CT, the physics of ultrasonic CT 

involves two steps; 

1. The measurement of the time-of-flight data for the wave passing through a number 
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of different paths through a test object, followed by 

2. The reconstruction of slowness coefficient values at all interior points from the 
data collected in step (1). 

Thus, acoustical tomography can be defined as a technique to determine values of a 
spatially varying acoustic parameter of interest across a slice through an object. Waves 
travel through the sample (from a series of source locations) to recievers at appropriate 
positions around the object. The source and the reciever locations are chosen such 
that the rays pass through as large a fraction of the object volume as possible. The 
relative positioning of the source and reciever should also conform to any requirement 
for regular positioning in the reconstruction procedure. As the source and the detector are 
moved around the object, a number of individual projections (or measured data sets) are 
obtained. Each projection is composed of a transmitted intensity pattern which depends 
on the path concerned. After completing a scan consisting of a number of projections, 
the distribution of the slowness coefficients within the object is reconstructed by one of 
the several reconstruction techniques. 

The major problem in applying UCT is that ultrasound does not travel strictly in lin- 
ear fashion and suffers scattering, ray bending and the diffraction at the interface of any 
dissimiliar material. In medical applications, however, UCT can be used safely in soft tis- 
sues imaging w*here deviation from linearity in propagation is minimal. Among industrial 
applications, practical material charecterisation problem such as study of residual stress 
distribution has also been tackled by UCT. With proper guideline for collecting projec- 
tion data and proper modification in the reconstruction technique, accurate quantitative 
evaluation of the material properties of the object slice can be accomplished. 
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1.8 Literature Review 


Nondestructive testing and evaluation (NDE) of composite materials posses a challange 
to both researchers and applied technologist. Ultrasonic methods, however, are most 
common and practice widely because of the low cost, relative ease in their use and the 
amount of information that can be obtained from them. In this section, a state-of- 
the-art review of the ultrasonic methods for NDE of composite materials is presented. 
This survey covers experimental techniques in nondestructive application of ultrasonics 
to characterize them and detect different flaw/damage found in them. 

1.8.1 Experimental Techniques 

A comprehensive review article on ultrasonic NDE of advanced composite was presented 
by Henneke [1990]. The article covers a wide range of works, applied to composite 
materials, which include different modes of wave propagation, material characterization, 
attenuation measurement, C-scan techniques, ultrasonic spectroscopy. 

The most basic application of ultrasonics is in characterizing any material by mea- 
suring the velocity of wave in it. Such well known techniques of calculating wave speed 
by measuring time-of-flight, are described in Bray and Stanley [1989] and Krautkramer 
and Krautkramer [1990]. A wave pulse is established at one surface of the material by a 
transmitting transducer. A second receiving transducer is located at the far side of the 
material to receive the signal after travelling through one thickness of the specimen. This 
mode of operation is called pitch-catch method. In the other method, called pulse-echo 
technique, the same transmitting transducer may be shutdown momentarily and the sig- 
nal is received after travelling through twice the thickness of the specimen. In both the 
cases time-of-flight is measured and subsequently the wave speed is computed. 

Wave speeds, in different modes, can be directly correlated with the material prop- 
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erties and vice-versa ( Kolsky [1963]). In isotropic materials, two independent velocity 
measurements are enough to determine the complete elastic moduli. The attenuation 
of any ultrasonic wave propogating through a medium gives in the information about 
the state in quantity of the medium. Attenuation is one of the main characteristics of a 
propogating wave that is highly dependent on the material and its properties. 

Tan chert and Hsu [1973] studied attenuation measurement of glass/fibre reinforced 
unidirectional and cross piled composite laminates under cyclically applied tension and 
compression loads. The attenuation increased with time but did not get back to its initial 
loading upon unloading. Attenuation was noted to be frequency dependent, increasing 
with higher frequencies and vice-versa. 

Store and Clarke [1975] observed that the main factors influencing ultrasonic attenu- 
ation were delamination, voids and state of the cure of the resin. Fibre volume fraction 
did not influence attenuation much. They performed C-scan on fabricated panels and 
attenuation was measured by comparing the originals with and without the specimen in 
the water bath. The attenuation was considered as the total loss and back faces end that 
associated with the propagation through the specimen. 

Hay ford et al [1977] and Meron et al [1977] established a correlation between ultrasonic 
attenuation and failure strength in composite specimen. 

Meron et al [1977] concluded that ultrasonic attenuation is a more sensitive indicator 
of the degradation of the properties due to exposure of the specimen in the hot water. 

Hayford and Henneke [1979] presented the results of a model to find attenuation 
changes observed (in real time) in tension test of graphite/epoxy specimens. They con- 
cluded that more transverse cracks do not occur in piles if the crack density with some 
minimum saturation spacing is readied. The changes in attenuation depends on the 
spacing of the transverse cracks which acts as a diffraction grating for the passage of the 
ultrasonic beam. 
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Williams and Doll [1980] measured the ultrasonic attenuation in graphite/epoxy spec- 
imens that were subjected to compression-compression fatigue. A through-transmission 
mode was adapted and, for measuring attenuation, the amplitude of the recieved pulse 
was compared with that of the input pulse at four narrow band centre frequencies of 
0.5, 1.0, 1.5 and 2.0 MHz. They also noted that attenuation to be frequency dependent. 
Group velocity, however, was found to be frequency independent. 

Srivastava and Prakash [1987] meaured attenuation in glass fibre zirconia hybrid com- 
posite by through-transmission immersion technique using a pair of 10 MHz probes. At- 
tenuation was measured by comparing the amplitude of transmitted probe with that of 
the incident probe. The result showed a bilinear relationship to be existent between 
pre- fatigue ultrasonic attenuation and the fatigue life. 

As defined earlier, a C-scan is the presentation of variation of the ultrasonic atten- 
uation of features(s) as a function of the transducer in two dimensions relative to the 
surface of the examined specimen. More recent developments in C-scan the various 
imaging techniques to improve data presentation. Multiple thresholds are set and either® 
various shades of gray or different colours are assigned to the different threshold levels. 

Hagemaier et al [1970,1971] inspected boron, glass and graphite reinforced polymer 
matrix composites using ultrasonic C-scan. They observed that through transmission 
techniques had the advantage of not being affected by surface roughness, surface contour 
etc. They also noted penetration at higher frequencies to be inadequate. However, 1.0, 
2.25 and 5.0 MHz were suitable for through transmission. They found that the C-scan 
was quite effective for locating planner defects and suggested that refrence specimen with 
built in discontinuities were essential to perform correct inspection. 

Mod and Stephenson [1971] inspected boron/epoxy with known defects, using through- 
transmission ultrasonic C-scan with infocused probes. They found it neccessary to per- 
form scans with different amplification and threshold levies to display all artifacts. They 
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observed that complete loss of a through-transmission original was typical of a delami- 
nation or air void. With higher setting of the threshold, they observed that even normal 
signalhad considerable fluctuations. Ultrasonic C-scan has also been used to study dam- 
age initiator and development in composites (Sturgeon [1977], Nahas [1985]). 

In recent years attenuation has been given to the development of automated C-scan 
technique. With the advent of composite materials, requirements have emerged for the 
development of the new inspection procedures. Mahoon et al [1982] discussed how the 
necessity of a skilled personnel and the variability of the personal judgement was overcome 
by automatic scanning system developed by British Aerospace. The scanning system has 
facility of uianipul.it ing probes in an water-squirt er system that controlled data collection 
and analysis. 

C-seans of tubular graphite epoxy specimen used in aerospace industries were per- 
formed by Rogovsky [1985] . In this investigation, flaws were simulated in different layers 
of tube and a method called “multiple reflections” was suggested where the number and 
amplitude of multiple reflection (from the boundary of the bonds) were considered. He 
concluded that an optimal selection of the sequence and number of signals for the flaw 
detectors gating yielded good detectability of the simulated flaws. 

Preuss and Clark [1988] used time-of-flight (TOF) ultrasonic C-scan for detection, 
sizing and characterization of defects in carbon 'fibre composite components. They ob- 
tained the depth wise information by recording the position of the defect echoes in time 
on the A-Scan. The C-Scan was constructed by using this TOF information rather than 
the back wall echo amplitude. They, however, applied the technique on 56-ply thick 
carbon/fibre composite panel with damage from impact and fatigue loading. 

Miller [1990] summarised a quantitative ultrasonic nondestructive technique applied 
to evaluate impact damage and porosity in graphite/epoxy composite laminates. Quan- 
titative imaged based on the slope of the attenuation coefficient measured as the function 
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of frequency over a broad bandwidth were obtained using a phase insensitive acousto- 
electric recieving transducer. Low velocity impact yielded increased values for the slope 
of the attenuation. 

1.8.2 Acoustic Tomography 

Computed tomography (CT) in particular X-ray CT, is a powerful NDE tool that has 
seen a rapid development in the past two decades. CT uses Computer reconstruction 
to combine the information obtained from multiple projections. The quantitative CT 
information is displayed as a reconstructed slice plane of ther part. Thus, in addition to 
qualitative evaluation, the CT data values can provide accurate quantitative measure- 
ment of material characteristics (density) in small volume elements of the object. 

Although CT is relatively new, most of its mathematical basis traces its origin to the 
pioneering work of Radon [1917] and his work on the reconstruction of funtions from the 
line integral data in the early part of the century. He showed that internal information of 
a three dimensional object (represented by a mathematical function) could be determined 
from the values of a set of line integrals of the mathematical functions. This work was 
expanded to radiography by Cormack [1963]. His work included both the mathematical 
approach and experimental verification. However, it remained until the invention of CT 
scanner by Hounsfield [1973], for the technique to be practically useful. 

Kak [1979] reviewed the major developments in imaging with computed tomography 
using X-ray, emission and ultrasound sources. The background and theory of the CT 
systems, types and applications of the CT systems can be found in the book by Herman 
[1980]. 

D.PJ arisen et al (1994] used Lamb wave tomography for finding out damages in 
advanced composite laminates. They studied two thin film laminate samples consisted of 
16 layers of carbon-fibre reinforced epoxy with a quasi isotropic lay-up configuration and 
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had been previously loaded in a biaxial test method until failure. The resulting damage, 
in the form of fibre failure, matrix cracking and delamination, has been imaged using a 
novel Lamb wave immersion tomography techniques. Images created with this method 
were correlated with images obtained with C-scan techniques. Both C-scan and Lamb 
wave tomography were able to identify clearly regions of damage in the two samples. 

William Wright et. al [1997] studied an entirely air-coupled inspection system using a 
pair of micro machined silicon capacitance transducers to image defects in the plates of 
different materials (0.7 mm to 2.22 mm thicness) using air-coupled Lamb wave tomog- 
raphy. Filter back projection algorithm was used in a form of difference tomography to 
reconstruct images of defects upto 10 mm diameter machined in aluminium and perspex 
plates, as well as in samples of carbon fibre reinforced polymer (CFRJP). The technique 
was able to resolve non-central defects as well as multiple flaws within the scan area. 

D.P. Jansen and I).A. Hutchins [1991] utilized ultrasonic Rayleigh and Lamb waves 
to reconstruct tomographic images of defects within metallic samples they used through 
immersion techniques, Rayleigh waves were used for the identification of surface struc- 
tures on thick samples, while Lamb waves are used to detect structures in thin sheets, 
they obtained tomographic images which gives information concerning the location and 
nat ure of the defects. 

D.P. Jansen et al [1993] obtained ultrasonic tomographic images from objects of rect- 
angular cross section by scanning spring-loaded transducers across the sample phases. 
The complete system was under the control of the microcomputer, resulting in a en- 
tirely automatic data collection system. Selected images obtained from this system are 
presented for metal and rock samples, as well as samples with artificial anomalies. 

Most recently Datta [1995] had developed a methodology in ultrasonic NDE for iden- 
tification and reconstruction of defects in fibre composites. He used Algebraic Recon- 
struction Technique (ART) for reconstruction. 
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1.9 Thesis Layout 


* Chapter 2 gives the detail of the data collection and the algorithm used for 

reconstruction. 

* Chapter 3 gives the detail of the experimental ultrasonic set-up used for scanning 
the specimen using ultrasonic CT scanner. 

* Chapter 4 describes the details of the program used for reconstruction, simulation 
and color coding. It also gives the reconstructed images for simulated data for validation 
of CRP and discussions on the results of the simulated images obtained from the simulated 
data. 

* Chapter 5 provides discussions on the results based on the real data. 

* Chapter 6 summarizes the conclusion and suggest scope for the further works. 
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Chapter 2 


Theoretical Formulation 


In this Chapter, the CT related computational and mathematical procedures underlying 
the data collection, image reconstruction, display are discussed. 

2.1 Preliminaries 

J'»(r,<p)dl=p(s,6) (2.1) 

where p[s,0) is called the projection data for the tomographic algorithm, and it is the 
line integral of function along the line specified by s and 6. (see Fig 2.1) 

Now the aim is to recover the value of p (r, <f>) ,the slowness coefficient, if a set of several 
p- values measured along various chords is given. This is the fundamental problem of 
CT and CBP has been used in the present study for that purpose. The p values can be 
suitably normalised to get the density values, if so desired. 
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2.2 Data Collection Mode 


The image processing methodology requires the slowness data to be collected by an array 
of ultrasonic recievers for the reconstruction of the function fi (r, <£). 

In this study the mode of collection is parallel-beam (PBG) mode (Fig.2.1).This 
system consists of several pairs of ultrasonic transmitters (S) and ultrasonic recievers 
(D) which can scan the object completely. 

The S-D pairs are spaced uniformly and the object can be rotated to get the data for 
different views. The line SD represents the path of the data ray or the chord along which 
V values is found out. 

The perpendicular distance from the center of the object to the path of the ray is 
denoted by s. Several SI) pairs could be used to get the p data. The object table is 
rotated to get the several sets of P for different values of 6. 

2.3 The CBP Algorithm 

The convolution back projection (CBP) algorithm has been described in detail by Herman 
(1980), In this section we review CBP briefly as reported earlier by Munshi (1992). 

Fig. 2.1 shows the data collection geometry for a parallel beam ultrasonic scanner. 
The object function /(r, 4>), is represented by a unit circle and on (of many) data rays is 
lrpiehcnled by SD. The ray indices are s and 6, where s is the perpendicular distance of 
the ray from the object center, and 0 is the angle of source position (or object rotation). 
The CT data is denoted by p(s, 6) given by 

p(M) = [ f(r,<fi)dz (2.2) 

J SD 

Here, z is the variable of integration along the chord SD. The ultrasonic setup collects 
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the projection data p(s. 6) for many values of s and 0. 

The ‘‘Projection- Slice” theorem [Herman (1980)] states the equivalence of the two- 
dimensional Fourier transform of / (r, (j>) and the 1-dimensional Fourier transform of 

7 >(.s\ 6) with respect to s. Symbolically 

p(R,6) = f(Rcos9, Rsin6) (2.3) 

Where the symbol A represents the Fourier transform and R is the Fourier frequency. 
A two-dimensional Fourier inversion of Eq.(2.4) leads to the well known tomographic 

inversion formula. 


r tt /*4‘OC' 

= / / p(R,0)e i27,Rrcos ^- o) \R.\dRde (2.4) 

J o J- oc 

Where, 

p(R,0) = J p{s,e)e-' 2 * Rs ds. (2.5) 

We note* that the inner integral in Eq.(2.5) is divergent. A practical implementation 
of Eq.{2.5) incojporate.s the replacement of the factor |/f| by \R\ H'(R), where W(R) is 
a suitable window function which vanishes outside the interval [-J? 0 R c ] . Here |f? c | is 
the Fourier cut-off frequency. Normally IV(R) is an even function of R . Thus Eq.(2.5) 
takes the approximate form, 

J /-7T #*4"00 

' / p(R,8)e i3nRrco,ie -^ |i?| W(R)dRdB (2.6) 

0 J *~oo 

If p(RJ) also vanishes for [i?| > R c , then the reconstructed f, denoted by /, agrees 
exactly with / the following window, 
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W(R) = { 


1 , \R\ < R c 
0 , |jR| > R c 


(2.7) 


An implementation of Eq.(2.7) in the spatial domain exploits the convolution theorem 
of Fourier transform so that the reconstruction, / , is given by, 


Where. 


and. 


r 7T r + OO 

f(r,<f>)& / p(s,6)q(s' — s)dsd8 

Jo J- oc 


/ +oo 

\R\ W(R)e ,2nR3 dR 

•00 


( 2 . 8 ) 


(2.9) 


s' = rco${6 - 0) (2.10) 

The index s' is of the data ray passing through (r, <&), the point being reconstructed.The 
inner integral in Eq.(2.8) is a 1 -dimensional convolution, and the outer integral (corre- 
sponding to the averaging operation over 6 ) is termed as back projection, and hence, the 
name convolution back projection for this perticular implementation. The CBP method 
is also known as the filtered backprojection algorithm due to the ’filtering 7 of the Fourier 
transform of the projection data, p, by the window ( or filter) W {R) in the initial stages 
of formulation as given by Eq. (2.6). The function q(s) known as the convolving function, 
is evaluated once and stored for repeated use for different views ( or different angle 6 ). 
So, for a given point (r, <f>), the inherent error JSj, in the CBP implementation, Eq.(2.8) 
is. 


Ei(r,d>) = f(r,<j>) - f(r^) 


( 2 . 11 ) 
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Where / and / are given by Eq.(2.4) and Eq.(2.6) respectively. This error is strictly 
due to infinite cut-off, R c , of the Fourier frequency and is precisely zero if the projection 
data happens to be band-limited and the cut-off frequency is chosen to be the highest 
frequency contained in p. In general, to avoid aliasing artifacts [Herman (1 980)], it is rec- 
ommended to choose, 


R c = 1/(2 As) 

Where As is the spacing of the data rays. 

It has been shown (Munshi(1992)j that E\ at a given point (r, <p) in the object cross- 
section. is given by 


£iM)w*(H'"(0))(V 2 /M)) 


( 2 . 12 ) 


Where, 


H’"(0) 


dHV(R) 
c m 2 


R = 0 


(2.13) 


V 2 / is the Laplaciari of /, 

and. k is a constant depending on the data-ray spacing. Eq.(2.12) is valid for objects 
having certain smoothness properties provided the data is perfect as per Eq.(2.2). 

The error £j represents the point-wise theoritical error in reconstruction, and it also 
obvious that the I.aplnnan of f(r,4>) has to exist for the predictions of the theorem to be 
valid. For points in the cross section, where V 2 / does not exist, the linearity between E l 
and W’"(0) is distured. If V 2 / is zero, other errors will be more dominant. The Laplacian 
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is zero for smooth regions of the objects, while it does not exist for rough edges. For 
simulated objects, E\ can be calculated from the original image, but for real objects the 
distribution is unknown, hence the error in reconstrustion cannot be claculated directly. 
This fact motivates and indirect representation of error. It has been reported earlier 
(Munshi et al [1997], Davis et al [1996]) that, for a given data set, sharpness can be used 
as an indicator of the behaviour of error, arising due to the choice of the filter function. 

If the image consists of a single point then the sharpness parameter is defined as the 
values of the reconstructed grey level at that point. For a general image, the sharpness 
parameter corresponds to, N MAX , the maximum grey level (linear slowness coefficient) 
in reconstruction. 

In the present work l/NMAX (inverse of sharpness) along with Eq.(2.l2) has been 
used to represent the physical nature of the object cross section. 
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Chapter 3 


Experimental Set-up and 
Procedure 

3.1 Present Set-up 

This Chapter describes the details of the experimental set-up, used for this study, to 
perform automated ultrasonic tomography. This involves generation of the ultrasonic 
pulse, conversion of recieved signal into digitized form in the real time and the storing of 
the data, 'this set-up gives the Tomographic images of the specimens. 

The present automated ultrasonic system, consist of a scanning tank, a high resolu- 
tion type Ultrasonic Flaw Detector (UFD), transducer, a high speed analog to digital 
converter (A-D) board and other accessories for precision movement of the transducers in 
two dimensions and data collection. The software controls the precision movements of the 
probes using stepper motors, digitize the recieved waveform and store it on the hard disk 
of the computer. The A-D board and PC-AT 386 machine, can digitize the ultrasonic 

signal in real time. A brief description about the set-up and the experimental procedure 
% 


23 



for automated probe movement and data acquisition is outlined in the following sections. 


3.2 Components of Experimental Set-up 

(1) Scanning Tank 

The purpose of the scanning tank is to house the stepper motors, the probes and the 
composite specimen. As is the common practice, the probes and the specimens are 
submtTjV'd in water which acts as a coupling medium. The side walls of the coupling 
tank are made of transparent perspex. Two lead screws, each having a pitch of 4mm, 
are fixed in mutually perpendicular directions on a frame located at the top of the tank. 
Each lead screw, mounted on frictionless ball bearings, is independently run by a seperate 
stepper motor. The probe holder is mounted on a mechanical base acting as a nut of the 
lead screw'. The? schama tie diagram and a photographic view of the set-up are shown in 
the Fig. 3.1 and Fig. 3.2 respectively. The probes can be moved very accurately as the 
lead screws are of zero backlash grade. 


(2) Ultrasonic Flaw Detector 

For generating ultrasonic waves and collecting the recieved signal an ultrasonic flaw 
detector (Krautkramer-Branson make, USIP 12) is used. The generated wave by the UFD 
can either be used in narrow band region or broad band region depending on the type of 
application. The UFD can be used with a single probe for pulse echo method or with two 
probes for through-transmission method. The recieved signal can be obtained from the 
UFD in rectified or in radio frequency (RF) mode. The gain of the recieved signal can 
be adjusted by fine controls to the desired level to reduce noise and avoid any saturation. 
The time base can be adjusted to clearly examine the successive echoes. Apart from these, 
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Figure 3.1: Sehamatic Diagram of Ultrasonic Setup 
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Ultrasonic flaw detector 



! I Motor Power Supply j j 1 

1 iP"~fr J 11 

Figure 3.2: Photographic view of the Ultrasonic Setup 


the UFO is also furnished with facilities such as noise reject, selection of two simultaneous 
time gates with threshold values, alarm facility and distance amplitude correction (DAC). 
The probes are connected to the UFD through standard LEMO connectors and the RF 
signal can be tapped from the UFD through a standard BNC socket for digitization. 


(3) Computer 

The set-up is run through a mierooomputei of PC- AT 386 SX family which is interfaced 
with the hardware i.e., stepper motor driver board and ultrasonic flaw detector via A- 
D converter. The scanning procedure is completely controlled by the computer when 
the scanning code is executed. The computer, with the storage capacity of 300 MB, is 
also networked with the central super mini computers so that the collected data can be 
directly sent to the systems for post processing. 
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(4) Stepper Motors 


The stepper motors used in the present set-up are of type STM 601 with 12 v input. The 
torque capacity of the motor is 2 kg-cm. The maximum step rate is 1023. The motors 
can he moved either in clockwise or counterclockwise. The motors are operated by the 
pulses sent from the driver board, the number of clock pulses for one revolution being 
200. Poi one revolution of the screw the probes move by 4 mm. 


(5) Stepper Motor Controller 

The POL-211 stepper motor control card, (manufactured by Dynalog Micro Systems), 
is used to drive the stopper motors in precise increments, direction, and speeds. Very 
accurate repeatable movements of the motor can be achieved. With this card a maximum 
of three motors can be controlled simultaneously. Each independent stepper channel 
consists of a 8039 processor that can execute a set of motion control commands. Every 
8039 processor also comes with an on board ROM where the control program for the 
processoi resides. All commands set by the PC are executed by the processor running 
the relevant control program. Once the command is set, the PC is set free, thus optimum 
uilisation of PC resources being possible. 

(6) Analog to Digital Converter 

The recieved RF signals from the UFD is digitized using a Matec A-D board. The board 
is capable of digitizing the RF signals with a maximum sampling rate of 100 MHz and 
maximum data acquisition flexibility through a 64k on-board high speed memory buffer. 
The various board functions, under software control, are input channel selection, AC/DC 
coupling, input voltage range, RF/Video mode, sampling rates, trigger selection, clock 
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control threshold phase and level, board selection and interrupt e. 5 

high speed data equisition memory is mapped directly into PC memoi v 

wave form has been captured, the PC can transfer data off the board at 

bit) or 3 MHz rate (16 bit) using simple memory block instructions. This b 

freedom from the static architecture of a stand alone instrument and the con, 

D100A oscilloscope is used to examine the waveform for selecting the appropi 
parameters such as delay, gate start, gate length etc. 

(7) Transducers 

Tlif transducers (probes) used in the present investigation is E16B which generates broad 
band longitudinal waves. The diameter of the transducer is 10mm. They are connected 
to the IT'D by water proof connectors. The fixtures hold the transducer in vertical 
position. 

3.3 Specimen Preparation 

The present investigations were performed on composite with polystyrene as the rein- 
forced material and perspex as the matrix material. The specimen are of circular shapes 
with diameter 60 mm and following are the description of the specimen used. 

SPECIMEN 1 (El): 

A 40 mm circular inclusion at the geometric center is made. First 40 mm hole is 
made in the perspex and the paste of polystyrene is filled inside the hole. The specimen 
is heated inside the oven to harden the resin. This specimen is as shown in Fig. 3.3. 


SPECIMEN 2 (E2): 
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In this specimen a 15 mm circular inclusion is made. The procedure of making this 
specimen is same as for the previous specimen. This specimen is as shown in the Fig. 3.4. 

SPECIMEN 3 (E3): 

In this snedmon a 25 mm x 25 mm square inclusion is made at the geometric cen- 
ter. The procedure for preparing the specimen is same as for the specimenl. The above 
specimen is as shown in the Fig. 3.5. 

SPECIMEN 4 (E4): 

4) This specimen is madeup of two circular inclusion of different sizes. One inclusion 
is of 15 mm diameter and other inclusion is of 10 mm diameter and both these inclusions 
are at the diameter of the specimen. The spacing between these two is 20 mm. The 
procedure for preparation is same as that of specimenl. This specimen is as shown in 
the Fig. 3.6. 

These circular specimen are housed inside a square wedge as shown in Fig. 3.7. with 
a circular hole of 60 rnm diameter. This is done to ensure the normal incidence of the 
ultrasonic wave so that the transmission through the specimens is maximum and the data 
for as many views can be collected. Earlier when square wedges are used, projections 
were limited because for each projection a separate set of wedges were required. So this 
problem is over come by choosing circular specimens and a wedge as shown in Fig. 3,7. 
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Figure 3.5: Specimen 3 


Figure 3.6: Specimen 4 
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Figure 3.7: Wedge 

3.4 Experimental Procedure 

The specimen to be investigated is cleaned and put on the rotatary table inside the square 
wedge. The probe is brought by moving the motor to the starting position as shown in 
the Fig. 3.1. It is ensured that the probe and the specimen are completely immersed in 
water. The distance between the probe surface and the specimen is so adjusted that they 
are at minimum distance without touching each other. The gain adjustment controls are 
so adjusted that the peak of the signal does not exceed 80% of the maximum screen 
height of the UFD to avoid saturation. The RF signals from the UFD can be viewed on 
the digital oscilloscope provided by the software of the A-D board. The gate perameters 
on the signal to be digitized are choosen and corresponding changes are made in the 
data file. Once the initial settings are done scanning of the specimen can be started by 
executing the master scanning code. The data collected is stored in the files on the hard 
disk in append mode. The scanning can be interrupted at any intermediate stage and 
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retailing ( an ^ f * nf m f ^ u * point where it was interrupted. The peak amplitude of 
th< Mgnal in tin* ft** 1 ! 111 11 % ^°*nain are extracted during the scanning process. To reduce 
the tine 1 dependent rS|iltt> amplitude detection at any location is repeated over a 

number of and t!o a>;o value is taken. 
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Chapter 4 


Data Analysis 

4.1 Numerical Implementation of the Programme for 

Discrete Data 

In practical application of reconstructions, the measured data correspond to the estimates 
of p(s,0) fi*t numerous values of s and 6. And the aim is to find the function j u(r,«j>) . 

Further, for the simplification of the reconstruction algorithm the projections are 
sampled uniformly in both s and 0. Therefore, for each view there are M equi-spaced 
rays As apart, ami N such views A# apart. We define 

W .=(M-W } ForMI 

M~ = ~{M - l)/2 

/ 

M* = (A//2) \ 

> F or even M 

M~ = -(M/2) 

4 
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In order to ensure that the collection of rays spesified bv 


{(mAS’.nAtf) : M~ <m< M + , 1 < n < N} 

rovers the unit circle from a complete range of directions, we choose 

Aft 7 r/JV 
and As ... l/'A/ 1 

4.2 Programme for Simulating Data 

The simulated data has been generated by Kumar [1996] for the following cases: 

1 , Knit circle (SI ). 

2. Knit circle with a hole of 0.3 unit radius (S2). 

3 Annular hole in the unit circle (S3). 

4. ! mct.’.'ir cross-section with nine holes of different sizes at different locations 

(Sis 

The number of rays, Mi AY and the number of views, N ANGLE can be varied as per 
the need The pnyr.mm.e written is a general one. For circular geometries, the region 
has beet; divided into three sectors of radius r%, r 3 , r 3 and the absorption coefficient 
/i . /,*. /,i respectively. For rectangular geometry also the programme will work for any 

number of holes, 

4.3 Programme for Reconstruction 

The cumulation back projection algorithm uses the data of parallel beam geometry to 

reconstruct the image of a cross-section. 
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The project ion data has been kept in the input file after applying the correction due 
to wedges which are used so that the incident ray must be normal to the test specimen. 
If the incidence is not normal, large part will be reflected then the transmitted, so a 
very v-cak signal at the ultrasonic reciever. The programme writes the reconstructed CT 
numbers in the result file which besides CT numbers also contain maximum, minimum, 
cross sectional!}’ averaged and the centerline averaged LITF values. In the present work 
eight Allens have been used to study the tomographic image. Beside the Ramachandran- 
Lakhnunarayan filter , sine, cosine and five Hamming filters have been used for different 
values of B, be B 0,54,0.75,0.80,0.91 and 0.99 . Where Hamming filter with B = 0.54 
is called Hamming smooth filter and Hamming filter with B = 0.99 is called Hamming 
shuip filter. The following filters are given by : 


UV/u 


< 


1 

Ham - Lak 

tm{ nliJ-JH.) 

Sin 

rm(rrfi/2H t ) 

Cos 


B <4(1- B)am{7tRJR c ) 


Hamming 
(0.5 < B < 1.0) 


The details »( the filters tire listed in Table 4.1with their Fourier-Space origin second 

. I. :a, |H , "{0)|, 
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[ Filter code 

B 

W'//(0) 

Ram-Lak 

- 

0.000 

Sin 

- 

0.083 

Cos 

- 

h 0.250 

h99 

0.999 

0.001 

h91 

0.917 

0.083 

~y_~h80 

0.800 

0.200 

h75 

0.750 

0.250 

h54 

0.540" 

0.460 


Table 4.1: Filter details 


4.4 Programme for Displaying Reconstructed Image 

and Colour Levels 

Fnt the put pose of displaying the image, the CT numbers from the result file is read and 
the ;$ppD»pn!e rutom is assigned corresponding to that. At the position of each number 
a small rectangle is diawn and filled with gray level corresponding to that CT number. 
The image produce is a pixel Thus by generating pixels for each CT numbers at their 
jiuMtions the> ran he graphically diplayed. The colour code is printed adjacent to it for 
quick l eft Piter 


4.5 Colour Code 

The rob mi rode is must to understand and extract information from the colour images 
produced in the present work. The natural choice for the colour code was our own solar 
spectrum IBCYnli' in the scale 0 to 255. Number 0 is represented by ‘Voilet’ ann 
255 by ‘Red’ and the intermediate colours are linearly interpolated to represent other 

numbers. 

R. i tun turn of tomographic data results in a matrix in which each element corre- 
sponds to the magnitude, of the material property being observed, at that point. If we 
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arc able to view this matrix in the form of a coloured image, each colour in the image 
will ( i>:it‘>{ioi!(l to the different observed values of the material property, thus providing 
un with a view to the inside object under observation. The VIBGYOR pattern provides 
a good mode of visualizing this variation in the material properties of the object. This 
visualization had been extremely simple if we had a fixed number of values for the ob- 
sei vahie material property and we had a unique identifier for each colour in a VIBGYOR 
spectrum. I hus the problem can he broken up into two parts: 

1 Identifying the RGB values for each of the colours required along the VIBGYOR 
spectrum. 

» Assigning the elements of the matrix to one of these colours based on a linear 
st almg of all the values in the matrix. 

It is a common practice to represent a colour as a set of three numbers, which cor- 
respond to the intensities of Red, Green and Blue respectively, that when mixed would 
produce tire colour, To tackle the first problem , we took the RGB values of the seven 

main colours of the rainbow as: 

* Voile* : MM, 0,211 

* Indigo: 75, 0, 130 

* Blue ■ 0. 0, 255 

* (been . 0. 255, U 

» Yellow: 255, 255, 0 

* Orange: 255, 165, 0 

* Red : 255, 0, 0 


In each case (J corresponds to the minimum possible intensity and 255 corresponds to the 
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maximum possible intensity of the colour (Red, Green or Blue). The program prompts 
the user to enter the number of colours desired along the VIBGYOR spectrum. Keeping 
the values of the major colours as given above, the other values are computed by linearly 
intei pointing the Red. Green and Blue colour components between these major colours. 

If 50 colours are inquired to lie between Blue and Green, then the Red component of 
all these colours is assigned a value of zero and the Green component is calculated as 
tin era, - 0 d {(255 - 0) * if 50), and similiarly the Blue component of these colours is 
rah uiatrd as Blur t 255 -t ((0 — 255) * i/50). 

Hr** second nb pi* dd. in is also tackled by using linear interpolation. The minimum 
and maximum values otvming in the matrix are determined first. If the number of colours 
lequned by user is «, then n ~ 2 values lying evenly spaced between the two limiting 
value> ate computed Now we have n values corresponding to each of the n required 
colours. While displaying the image the colour the colour to be assigned to each pixel is 
determined by eompaiision with these ti values. This can be better understood in light of 
the following example Suppose the user asks for a 256 colour long VIBGYOR spectrum, 
|**l the smallest value found in the matrix is r mtri and the largest value is V maz . An array 
of 256 numbers is created and V mm and Y m a* axe placed in the first and the last place, 
j , . , These and the other 254 values can be computed using the formula: 

\\ , Y m n *f {Km* - Knin) * t/256 0 < i < 255 

4.6 Validation of CBP from Simulated Data 

The simulated data of Kumar (1990) have been studied for four different samples as 
mentioned in the sertum (4.2). The results have been summarised in table 4.2 to table 
4 r, The original and the reconstructed images are shown in Fig. 4.3 to Fig. 4.10. Re- 
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construction has been done for h54 and h99 filters. Unit circle is easiest to reconstruct 
and gives the most accurate reconstruction. These samples has been chosen to represent 
tin* flaws in the test specimens with increasing complexity. 

By looking at the reconstructed images, we conclude that the results obtained are 
good. Tim quality of reconstruction of simulated objects are very good because of ab- 
sence of am error but the error in GBP algorithm. The gray level is symmetrical for 
symmetrical objects. 

The gray -level do not changes sharply at the edges as in the actual objects. Thus the 
edge effects are visible. But the density of the reconstructed section of simulated data is 
unifor m at 4 to ft pixel away from the edge of the pixel boundary. 

T he above analysis on the simulated objects shows that CBP algorithm can give good 
results for the symmetric as well as for noil-symmetric objects. 

Pig -1.1 ami Fig 4.2 shows the comparative performance (1/NMAX) of filters for 
simulated samples Kefering to section 2.3, for a given cross-section distribution 
the point wise error in reconstruction is proportional to the magnitude of kP/r(0), provided 
that the collected data are “perfect”, and the specimen characteristics around the point 
in question is such that Faj.(2.12) is valid [Davis et al. (1995)). So processing of a given 
data set by several filters results in a linear relationship between W7/(0) and Ei. The 
plots have been shown in from Figs. 4.1 and Fig. 4.2 . The plot for simulated specimens 
Si, 82 and S3 shows wide variation from linearity. Since these samples are smooth, V 2 / 
vanishes for the major portion of the region, hence, Eq.(2.12) no longer represents total 
error in this rase. Ftg. 4.2 show good linear fit because for this sample V 2 / exist 
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Kilt eis 

NMIN 

NMAX 

1/NMAX 

hr,.i 

25 

75 

0.013333 

h7o 

25 

76 

0.013157 

" hBO 

25 

78 

0.012987" 

“ m "~ 

25 

81 

0.012345 

h«)D 

25 

82 

0.012048 


Table *1.2: Summary of results for simulated specimen 1 


1' liters j 

NMIN 

NMAX 

1/NMAX 

h">4 ’ 

22 

75 

0.013333 

| h7f> 

'23 

76 1 

0.013157 

bSO 

22 

78 

0.012827 

1 1,;>l 

IS 

81 

0.012345 

[ h*K» 

' ' 

83 

0.012048 


t able >1.3' Summary of results for simulated specimen 2 


j Filters [ 3 

s’MIN" 

"NMAX” 

1/NMAX 

ha 1 

24 

77 ] 

0.01298 

, i 

hr:. . 

24 ‘ 

76 

0.01315 

i hsu ‘ 

24 

7<5 ” 

0.01282 

♦ t 

. h9 i ; 

' 23 


”0.01219 

■ h'M i 

21 

81 

0.01190 


Table 4.4: Summary of results for simulated specimen 3 


fllteli 

NMIN' 

NMAX 

51 

1/NMAX 

0.01851 

hr>4 

•><) 

" 1575 “' 

21 

55 

0.01818 

KST“ 

"21 

56 

0.01785 

, — - 

20 

56 1 

0.01785 

hou 

r ““l9 

L tl 

57 

0.01754 


Table 4.5: Summary of results for simulated specimen 4 
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Figure *1.1: 1 NMAX plot for simulated specimens 



Figure 4.2: 1/ NMAX plot for simulated specimens 
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Fig. 4, 1 ; Original Image for the simulated specimen 1 


Fig. 4.2 : Reconstructed image with h54 and h99 filters 


Mm * 1 .0000 Max = 256.0000 

LAvg « 0.9727 AAvg * 0.8068 

File : h99 J, NRay : 37, NView : 60 
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Fig. 4.3: Original Image for the Simulated Specimen 2 






Fig, 4,5: Original Image for the Simulated Specimen 3 
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He. 4.6: Reconstructed Image with h54 and h99 filters 





Fig. 4.7: Original Image for the Simulated Specimen 4 
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Chapter 5 


Experimental results and 
Discussion 


Jj* the. ( |. 'j.'i’i *i ir troth-. obtained from the actual experimental data have been dls- 

oiv,-.*, 


5.1 Data l ' Keel 

1 fa 4 ,»f.t u *.»• ukrn with th** help of att ultrasonic setup. The details are briefly 

Mill. ill. s! }'.« 4 h*'le 

4 I V u^.tv :.!> u.u.Mni»m and mdevers are fixed horizontally at a distance of 9.35 
. m 1 V t tthidi is used is of 5 MHz frequency. This set of transmitter 
;.-»* > 4 <»- Jji-srd it! 4 hwmmtai direction at a regular step of 1 mm upto 90 mm. 
h : ? },. p, . *.v< ti vd on the rotating table which can be rotated with the help 

of .* M, }) p, , ^ ,.i tMmwlh Hen* manual rotation was used. The specimens are 

„, U n.4 in •» <!► gp-rs np'u l(Mi . n > • '< .m die specimen completely. 
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O The tint a is rollwlcd for four different specimen in which two axe radially symmet- 
ric. one is quadrant symmetric and the last specimen is unsymmetric. All the specimen 
ar«* made of perspex and the inclusion material is polystyrene. The acoustic impedance of 
t he>e t w o material is differ by 1051 ,The information of the different specimen is as follows: 


(i) SPECIMEN 1 (El) 

This specimen is a circle with radius 30 mm and a circular incision at the center of 
radius 20 mm The geometry of this specimen is shown in Pig. 5.3. This specimen is 
scanned » it h 90 rays and for 90 views. The distance between two consective rays is 1 mm. 

in) SPECIMEN 2 (E2) 

This specimen is also a rircle as the previous specimen but here the circular inclusion 
is of radius 7 5 mm The geometry of this specimen is shown in Fig. 5.5. Here also the 
spmmrn is scanned with 90 rays and for 90 views. Here also the distance between two 
ronsrrutisT rays is 1 mm 

(m) SPECIMEN 3 (E3) 

This specimen is circle with radius 30 mm and a square inclusion at the center of 
length 25x23 mm The geometry of this specimen is shown in Pig. 5.7. This specimen is 
scanned with 90 rays and 20 views. Here also the distance between the two consecutive 
rays is same as the previous specimen. 

(iv) SPECIMEN 4 (E4) 

This specimen Is also a circle with radius 30 mm and two circular inclusion at the 
diameter of Iff mm diameter and 20 mm diameter respectively and the distance between 
them is 23 mm The geometry of this specimen is shown in Fig. 5.9.This specimen is 
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alM> sramifd with 90 rays and 20 views. 


5.2 Discussion of the Results 

Thf i ( ,k -ii!ts of tin* experiments done on four different specimens with eight different class 
of fillers are obtained But we take only two filters Hamming 54 (h54) and Hamming 99 
fuihW* for discussion. As all the results with different filters are more or less similiar to 
these ri'MjIts o: lies between these two . Discussion on the results of the experiments is 

as follow* 


It; M '1 t iV! X 1 (Kl>: 

Hie original ami the reconstructed image of the specimen are as shown in the Fig. 
5 a and Fig 5 1 Tin- first part of the Fig.5.4 is reconstructed with Hamming smooth 
filter » e ha I and the second part is reconstructed with Hamming sharp filter i.e. h.99. In 
both the moHsitucted images the inclusion is clearly seen at the center of the specimen 
but which is of small M/e as compared to the inclusion in the original image. The reason 
foi this is a * the ultrasonic ray traverse a medium whose acoustic impedance is different 
from the one m which it is travelling, it bends around the inclusion due to which time 
of flight im te.wM h and reduc es the size of the inclusion. There are two circle visible of 
diffident radii which is due to the following two reasons. 

jai Wobble m the rotating table of the ultrasonic setup and 

(hi The mismatch in the centres of the specimen, specimen holder and the rotating 

table 

There is some error also because of the water in between the specimen and the wedge. 

(it) Ml « i Ml X 2 (Eli): 

The original ami the reconstructed image of this specimen are as shown in Fig. 5.5 
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and Fig. 5.6. This specimen is also reconstructed with two different filters h54 and h99. 
Here the inclusion of smaller radii as compared to the previous specimen is very clearly 
seen. Same explanation above is also valid for this specimen also. 

(iii) SPECIMEN 3 (E4): 

The original and the reconstructed image of this are shown in Fig. 5.7 and Fig. 5.8. 
This specimen is also reconstructed with two different filters h54 and h99. Here some 
lines in the star formation are clearly visible which are due to the lesser views in scan- 
ning the specimen. The inclusion at the center is not a square which is due to the sharp 
edges at the square corners and due to which ray bending effect is more. So from the 
reconstructed image defect can be detected but with some truncated edges. 

(iv) SPECIMEN 4 (E4): 

The original and the reconstructed image of this specimen is shown in the Fig. 5.9 
and Fig. 5.10. This specimen is also reconstructed with two different filters Hamming 
54 and Hamming 99. With Hamming sharp (h99) filter high frequency components are 
clearly visible which are not in the Hamming smooth (h54) filter. Here many small 
inclusions in place of two bigger inclusions are visible. The exact size of the inclusion is 
not determined becuase when ultrasonic wave travers two inclusions the path it travel is 
not straight and becuase of which error are introduced in the reconstruction. 


5.3 Comparative Study of Filters 

All the tour specimen are examined with five different Hamming filters. The results 
have inirnm*— 1 to table 5.1 to table 5.4. As it has been discussed earlier that, tor real 
objects the distribution is unknown and the error in reconstruction cannot be calculated 
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dnetlly, It has been reported earlier that, for a given data set, sharpness can be used 
as an imhratior of the behaviour of error, arising due to the choice of the filter function. 
Ha sharpness graphs for the real objects are shown in Fig. 5.1 and Fig. 5.2. The graphs 
slews a linear relationship with different intercept and slopes. 

Tin' thoiem which was given by (Munshi [1992]) has been already verified for X-rays 
f Mniislu et al [1993]). (Davis et al [1996]) and for X-rays and gamma rays by (Wells and 
Munshi {1994]l, But for the first time an ultrasonic set of data is used to verify the above 

fttenfeiu, 
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Table 5.1. Summary of results for Specimen 1 


Ai‘ I." . 

NMIN 

NMAX 

1/NMAX 


■56(18 

107363 ' 

9.31 x 10 -7r ” 

1.75 I 

■7733 

*111917 

8.93 x 10^~ 

i- vi 

-S2U7 

114162 ' 

8.75 x 10^ 

h ( u : 

j 

-12*57 ' 

120226 ' 

8.31 xr T 

}.!N , 

- I (M3 1 " 

124398 

8.03 x 10" 6 


Table 5,2: Summary of results for specimen 2 
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Table 5.3: Summary of results for specimen 3 
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Table 5.4: Summary of results for specimen 4 
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Min * 1 .0000 Max * 256.0000 
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Fite : h99.fi!> N Ray : 90, NView : 90 
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ted Image with h54 and h 99 filters 
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Fig. 5.4; Reconstructed Image with h54 and h99 filters 



Fig. 5.5: Original Image for the Specimen 3 
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Fig. 5.6: Reconstructed Image with h54 and h99 filters 
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Mg. 5.7: Original Image for the Specimen 4 
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i-jg, 5,8; Reconstructed Image with h54 and h99 filters 


Chapter 6 


Conclusions and Suggestions for 
Ini in re Work 


Fit >m the study <»f the reconstructed images of the specimens, following conclusion can 

he dlaWU 

• Circular defects at the renter with 35 mm or more diameter can be easily detected 

and sized , 

* Square defec ts at the renter with size 25 mm x 25 mm or more can be detected but 
, A .»> t sizing cannot be done. In general defects do not have sharp edges as compared to 
square so defects which are similiar to square can be detected and sized. 

♦ Two or more defec ts of different sizes cannot be detected very easily but in such 
type of specimen it can be said that there are two different acoustic impedance zones. 

* The results can further improved if more no of projections with more number of 
i ay*, are taken and with smaller diameter ultrasonic probes. 
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Future Work 


* Data ran he collected for specimen with small size inclusions of different materials. 

# A suitable method should be adopted to remove the wobble in the rotating table of 

f hr ultras mie M*tup. 

# A Miitablr method should also be adopted to match the center of the rotating table, 

sprt nnrii huldn and the specimen. 

* Ha> tt at ntj; can he applied. 
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